Continuous THz lasing from dipolaritons 



m 

o 

Oh- 
< 



K. Kristinsson,! O. Kyriienko/'^ T. C. H. Liew,i and I. A. Shelykh^'^ 

^Division of Physics and Applied Physics, Nanyang Technological University 637371, Singapore 
^Science Institute, University of Iceland, Dunhagi-3, IS- 107, Reykjavik, Iceland 

(Dated: April 5, 2013) 

We propose a scheme of continuous tunable THz emission based on dipolaritons — mixtures of 
strongly interacting cavity photons and direct excitons, where the latter are coupled to indirect 
excitons via tunnelling. We investigate the property of multistability under continuous wave (CW) 
pumping, and the stability of the solutions. We estabUsh the conditions of parametric instability, 
giving rise to oscillations in density between the direct exciton and indirect modes under CW 
pumping. In this way we achieve continuous and tunable emission in the THz range, in a compact 
single-crystal device, which is expected to operate at high temperatures. We show that the emission 
frequency can be tuned in a certain range by varying an applied electric field and pumping conditions. 
Finally, we demonstrate the dynamic switching between different phases in our system, allowing 
rapid control of THz radiation. 
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Introduction. The development of methods for genera- 
tion of continuous radiation in the terahertz range (0.3-3 
THz) is currently an important physical challenge py]. 
Terahertz radiation has a wide area of application, rang- 
ing from astronomy [2[ and nondestructive spectroscopy 
[3[ to security and medical imaging. Various families of 
devices for THz radiation generation have been theoret- 
ically proposed and experimentally tested. They include 
emitters based on the solid state Gunn diode [j] and the 
free-electron laser [5[ , which allow for high output power 
of the radiation, while suffering from a bulky size. Mi- 
croscopic devices being developed for THz generation in- 
clude semiconductor structures excited with femtosecond 
laser pulses, where oscillations of charge density induce 
terahertz emission by classical dipoles [Ifl]. However, the 
tunability of this source and its output power are highly 
restricted. Another microscopic emitter, proposed the- 
oretically in 1971 [7|], generates coherent THz radiation 
through repeated intersubband transitions across stacked 
quantum wells. These quantum cascade lasers (QCLs), 
realized experimentally in 1994 jSj, achieve high efficiency 
and relatively high output power (up to 600 mW [9| ) , and 
a typical frequency of emission lying close to the upper 
bound of the THz range. A limiting factor of QCLs is 
the cryogenic temperature at which they operate. Thus, 
the lack of universal THz emitters with high power, ef- 
ficiency, small size, which are easily- tunable and operate 
at relatively high temperatures stimulates the search for 
sources based on other operating principles. 

Polaritonics — an area of physics which combines con- 
densed matter physics and quantum optics 12|, ll3| — 
offers new possibilities. For instance, a THz emitter 
based on transitions between upper and lower polariton 
branches was considered in Refs. 14l4l6|. Another pos- 
sible scheme exploits the THz range transition between 
2p and Is states of the exciton, the former being a dark 
state and the latter coupled to the cavity mode [17|. A 



recent study suggests a bosonic QCL, with multiple THz 
photon emission from transitions between energy levels of 
exciton-polaritons confined in a parabolic potential [18|. 

This Letter builds on the subject of a recent proposal 
to use a dipolariton system to generate pulses of THz ra- 
diation through oscillations in density between spatially 
direct and indirect excitons [19]. The dipolariton system 
consists of a semiconductor microcavity with a double 
quantum well (QW) embedded in the center [20|, |21| (see 
sketch in Fig. [IJa)). In a high quality cavity one can 
achieve strong light-matter coupling between cavity pho- 
ton mode (C) and direct exciton (DX) in the left QW 
(Fig. mb)). In the same time, tuning electron levels of 
the left quantum well (LQW) and right quantum well 
(RQW) into resonance, one can achieve the strong tun- 
nelling coupling between a direct exciton and a spatially 
indirect exciton (IX) formed by an electron in the LQW 
and a hole in the RQW |22|-|24|. These strong couplings 
between three initial modes lead to the appearance of new 
eigenmodes of the system, which represent three linear 
superpositions of bare cavity photon (C), direct exciton 
(DX) and indirect exciton (IX) modes. They are called 
the upper dipolariton (UP), the middle dipolariton (MP) 
and the lower dipolariton (LP). 

In this Letter we show that account of nonlinear ef- 
fects arising from exciton- exciton interactions can qual- 
itatively change the behavior of the system and allow 
achievement of stable continuous THz emission with tun- 
able properties. Nonlinear it ies are known to give rise to 
bistability |25l427j when a single mode is excited with a 
coherent pump slightly above resonance, and multista- 
bility [28|, l29| in configurations where additional states 
are available. Aside from the possibility of switching be- 
tween different stable states |30M33| , parametric instabil- 
ities can achieve periodic oscillations in particle densities 
[3J, l35|. Here, we first investigate the property of mul- 
tistability under continuous wave (CW) pumping, and 
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FIG. 1: (Color online) Sketch of the system, (a) Illustra- 
tion of the optical cavity formed by two distributed Bragg 
reflectors (DBR) with two quantum wells (LQW and RQW) 
placed inside. Labels identify the cavity mode (C), direct ex- 
citon (DX) and the indirect exciton (IX), which are coupled 
by Q and J coupling constants and form dipolariton modes. 
The red arrow shows the THz radiation emitted by an oscil- 
lating dipole formed by dipolaritons. (b) Band diagram of the 
double quantum well tilted by applied electric field. The left 
QW is coupled to the cavity mode. Electron energy levels are 
tuned to resonance, with hopping constant J between wells. 



the stability of the solutions. We establish the condi- 
tions of parametric instability, giving rise to oscillations 
in density between the DX and IX modes. In this way 
we achieve continuous and tunable emission in the THz 
range, in a compact single-crystal device, which is ex- 
pected to operate at high temperatures [36|. We show 
that the emission frequency can be tuned in a certain 
range by an applied electric field. Finally, we demon- 
strate the dynamic switching between different phases in 
our system, allowing rapid control of THz radiation. 

The model. The system is represented by two coupled 
QWs placed inside an optical microcavity (Fig. [1]). For 
a coherently pumped cavity mode and resonant coupling 
between QWs, separated by a thin barrier, the system 
hosts dipolaritons — strongly mixed modes formed by 
cavity photon, direct exciton and indirect exciton modes. 
We use a generic Hamiltonian of the coupled system: 

?i =huJcCi^CL + hujDxb^b + hujixc^c -\ {a^b + b^a) 



Vdi b^'c^'bc + P{t)d) + P{t)*a, 



(1) 



where a) , b^ and c^ are creation operators of cavity pho- 
tons, direct excitons, and indirect excitons, respectively. 
The first three terms in Eq. ([!]) represent the energy 
of the bare cavity {hwc)-, direct exciton {fvjODx) and in- 
direct exciton ihwix) modes. The following two terms 
describe the linear coupling between modes; the first de- 
noting the Rabi splitting between the cavity mode and 
the direct exciton, hVt^ and the second being the direct 
to indirect exciton tunneling rate hJ . The sixth, sev- 



enth and eighth terms introduce nonlinear interactions 
into the system. Vdd and Vjj are the interaction matrix 
elements between pairs of direct and indirect excitons, 
respectively. The inter-species scattering of a direct ex- 
citon with an indirect exciton is described by the matrix 
element Vdi- The two last terms correspond to coherent 
pumping of the cavity mode with intensity \P{t)\^. Un- 
der CW pumping the time dependent pumping rate can 
be written as P{t) = Poe~*^p^ where hwp is the energy 
of the pump and Pq is its constant in time amplitude. 

Equations of motion for the macroscopic order pa- 
rameters defined as expectation values of annihilation 
operators (a^) = Tr{pai}, i = C^DX^IX^ can be ob- 
tained using the Heisenberg equations of motion for op- 
erators a^b^c and applying the mean field approximation 
{aittj ■ ■ • GLk) ^ {ai){aj) ■ - - (cLk). Additionally we perform 
the change of variables a^ -^ e^^^^^a^. The resulting 
system of equations reads: 
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{Vjj\{c)\' + Vni\{b)\'){c), 



(2) 



(3) 



(4) 



where we introduced the lifetimes of the modes re = 5 
ps, tdx = 1 ns. Exphcit references to mode energies 
have been removed by defining relative energies 6q = 
^c — ^DX and Sj = ujix — ojdx • The pumping term now 
reads P{t) = e'^^^'^P{t)/h. Under CW pumping we have 
P{t) = Poe~*^^S where Pq = Po/h, and Ap = uop — ujc 
is the relative pumping frequency. . 

The conditions investigated in this Letter correspond 
to high occupation numbers for which the first order 
mean field approximation is applicable. To test this we 
also derived dynamic equations written for higher order 
mean field theory (see Supplementary material, section 
A). Numerical solution showed no change in the results, 
confirming the accuracy of Eqns. (J2j)-(j4]). 

In addition to dynamics we can study the steady-state 
properties of the system assuming a CW pump with tun- 
able energy and intensity. It was shown for the case of 
a pumped mode with nonlinearity present in the system 
that several solutions for the occupation numbers of the 
modes as a function of pumping intensity are possible 
37H4Q]. Moreover, in the case of two coupled modes 
the system can exhibit either stable solutions or para- 
metrically unstable solutions characterized by continuous 
oscillations of the mode occupation numbers [39|. The 
stationary solutions can be found using the ansatz for 
the modes (a^) = ?/^^e~*^p^ which corresponds to har- 



(a) 




(b) 


JB 1 '^^^B 


Occupation numbers (10^ particles) 


Mj 11 1 


3.0 


/ Ndx T=0.57 PS N,x 


jw if k 


2.5 


A L A A A A A aTa ft A ft A^ft A ft A ft fl 


4P // 1 


2.0 


'\ 




^ k IX 


1.5 






mF\l\ /xs 


1.0 




I I I 1- 


■^^1^^^^- 


0.5 


V V 


y y u u 


^^^^^^^^^^^^ 


0.0 


\/"V"l/'V-V-V'V'V'U'U'b 



2 3 4 5 6 7 
Mp (meV) 



502 503 504 505 506 
Tlnne(ps) 



FIG. 2: (Color online) (a) Phase diagram plotted for indi- 
rect exciton population numbers on the vertical axis and the 
pump energy on the horizontal axis. Enclosed regions in- 
dicate parametric instability (PI). Elsewhere the system is 
stable. Lighter colors indicate higher pumping strengths are 
needed to reach relevant population numbers. We can iden- 
tify the resonance with the MP mode and the UP mode as 
dark regions. Both modes shift higher in energy with higher 
population numbers. The white region in the upper left cor- 
ner indicates population numbers that can't be reached at 
this pumping energy. Pumping conditions in (b) are identi- 
fied with a red cross, (b) Time evolution of direct and indirect 
exciton populations in the region of parametric instability. At 
large times populations oscillate harmonically. 



monic oscillations with the frequency of the CW pump. 
Inserting this into Eqns. (|2])-(|11) and eliminating the non- 
interacting photonic mode tjjc^ we obtain a coupled pair 
of equations for the occupation numbers: 



(-^1 - i7i + VddI^dxI + VDI\^pIx\)^I^DX 




- -hJ^Pix + Pi =0, 


(5) 


{-E2 - i72 + VDili^Dxf + Vii\iJixf)lpIX 




- IhJ^DX = 0, 


(6) 



where £^^,7^,Pi, i = 1,2 are dressed energies, decay 
rates and pumping strength, with explicit expressions 
given in the Supplemental material, section B. Eqns. (|5j) 
and (|6]) represent a system of nonlinear equations for the 
IX and DX occupation numbers with multiple solutions 
in certain ranges of the pumping intensity. Their sta- 
bility analysis can be performed similarly to the method 
used in Ref. [39|, shown in the Supplementary mate- 
rial, section B; one adds a fluctuation to the mean- field 
solution and through substitution back into the dynamic 
equations one can determine if the fluctuation grows or is 
suppressed from the imaginary component of its energy. 

Results. We model the system with linear coupling 
parameters chosen as hJ = hfl = 6 meV [21]. The 
DX-DX interaction constant can be estimated as Vdd = 
QEi)0?^/ S |4l|, where the direct exciton Bohr radius and 
binding energy are a^ = 10 nm and Ei) = 8 meV, respec- 
tively. Here 6* = 100 fim^ is the laser excitation area. 
The IX-IX scattering constant was taken from Ref. [42|, 
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FIG. 3: (Color online) Multistability under CW pumping. 
Blue curves correspond to direct exciton occupation num- 
bers, and red curves to indirect exciton occupation numbers. 
Thin curves indicate stability of the population numbers, 
while thick curves indicate the presence of parametric insta- 
bility, (a) The pumping energy HAp = — 1 meV is slightly 
higher than LP energy. Here, the bist ability is induced by 
the blueshift of the LP mode and the DX number (left axis) 
is much larger than the IX number (right axis), reflecting 
the small IX fraction in the LP mode, (b) The pumping 
energy HAp = 4.5 meV is slightly higher than the MP en- 
ergy. Here, bist ability is induced by the blueshift of the MP 
mode and the IX number is of the same order of magnitude as 
the DX number, reflecting their almost equal fractions in the 
MP mode. This regime facilitates large oscillations in dipole 
moment under parametric instability, (c, d) The pumping 
energy HA = 8.5 meV is slightly higher than the UP energy. 
The main plot shows bistability coming from the blueshift of 
the MP mode. The pump intensity is much higher comparing 
to (a) and (b), as we pump at higher energy. The insets show 
mult ist able regions governed by the blueshift of the UP mode. 



with the QW separation taken as L = 12 nm (4 nm tun- 
nelling barrier). The derivation of the DX-IX interaction 
constant is shown in Supplementary material, section D. 
Relative energies of the modes were chosen as HSq = — 3 
meV and hSj = 1 meV, the latter corresponding to an 
applied electric field of magnitude F = 0.91Fo, where 
Fq = 12.5 kV/cm is the field strength at which the DX 
and IX modes are resonant [43. The numerically per- 
formed stability analysis of solutions of Eqns. (J5j) and (J6j) 
allows us to plot a phase diagram of the system shown in 
Fig. [21(a). For certain values of pumping energy HAp and 
pump intensity, connected to the number of indirect ex- 
citons Njx^ solutions that are not stable are possible. In 
these regions we can have parametrical instability (PI). 
Multistability curves plotted for pumping energies 
hAp = —1,4.5,8.5 meV are shown in Fig. [3](a), (b) 
and (c,d), respectively. Stable steady state solutions are 
shown with narrow, darker lines, and thicker lines in- 
dicate parametric instability. One can see that the DX 
populations are reminiscent of the conventional one mode 
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FIG. 4: (Color online) Frequency of oscillations of indirect ex- 
citon occupation numbers plotted as a function of the applied 
electric field, for cavity mode detuning 6n = —1 meV and 
pumping energy Ap = 4.5 meV. Inset shows anharmonic oscil- 
lations in indirect exciton numbers for electric field F — 1.2Fq. 



polariton system bistability [25|-[27|, with regions of in- 
stability now replaced with parametric instability. This 
implies that population numbers do not decay to stable 
steady state solutions, but oscillate continuously in time. 
Solving numerically Eqns. ©-(Ill) we show long-standing 
beats of IX and DX occupation numbers under CW pump 
(Fig. Efb)). Note, that this behavior of the system is only 
possible due to the presence of interactions, while with- 
out accounting for nonlinearities beats of IX- DX density 
quickly decay in time [19|. The pumping energy was 
fixed to Ap = 4.5 meV and the CW pump intensity 
was linearly turned on to |Pop = 1.1 • 10^^ s~^. For 
long turning on times (> 5 ps) the system remains stable 
on the lower branches of bistability, with Njx = 2500, 
and Ndx '^ Njx- For short turning on times (< 2 
ps) the system becomes parametrically unstable, result- 
ing at large times in harmonic oscillations of the popu- 
lation numbers. The system can also be switched from 
the stable behaviour to the parametrically unstable one 
by applying additionally to CW pumping a short laser 
pulse, as shown in the Supplementary material, section 
C. This represents the high degree of control of the pro- 
posed system. 

Furthermore, we study the behavior of the system sub- 
jected to the variation of the detuning Sj between IX and 
DX modes, keeping the energy and intensity of the pump 
constant. This change alters the frequency of particle 
numbers oscillations. The dependence of the frequency 
on the applied field is shown in Fig. 21 with frequency 
spanning from 1.5 to 2.2 THz. For the values of F past 
1.1 Fo oscillations become strongly anharmonic, shown in 
the inset in Fig. IH which means that the system demon- 
strates pronounced multi-mode emission. When the elec- 
tric field is decreased past the limits of Fig. HI the time 
to reach harmonic oscillations increases rapidly, limiting 
the usefulness of this range. 

Discussion. We have shown that exploiting parametric 



instability we can achieve harmonic oscillations in the in- 
direct exciton population under CW pumping conditions. 
An important property of spatially indirect excitons is a 
non-zero dipole moment directed in the growth direction 
of the microcavity. The single indirect exciton carries 
dipole moment equal to do = eL, where L is the QW 
separation, while the total dipole moment of the system 
can be calculated as D{t) = doNix{t). In the case of the 
oscillating population numbers, this can be written: 



D{t) = doN^x cos(cjt/2)^ 



(7) 



where u = 27ri^ is the angular frequency of the IX pop- 
ulation oscillations, and Njj^ is an amplitude of the IX 
occupation number modulation. Thus, oscillations of the 
IX density induce oscillations of the total dipole moment 
of the system with frequency lying in the terahertz range. 
This results in the emission of THz radiation by the ar- 
ray of classical dipoles formed from dipolaritons excited 
over an area in the microcavity plane. 

To estimate the power of emitted THz radiation, one 
can treat the system as a classical dipole antenna, which 
has a far field radiation intensity given by |44| : 



/o 



67reoc^ 



37reoc^ 



(8) 



where eo is the vacuum permittivity, and c is the speed of 
light. The radiation intensity is directionally dependent, 
as I {6) (X sin^(^), where is the angle of emission with 
respect to the microcavity growth axis. Choosing the 
applied field as F = 0.91Fo, the frequency of oscillations 
is 1.75 THz, and the maximum IX number is Njx = 
2.4 • 10^. The total emitted power is then /q ^ 14 nW. 

It is important to note that the intensity is quadratic 
in the indirect exciton number. This can be explained 
by the phenomena of superradiance [45|, l46| — coher- 
ence of the oscillations in the quantum system causes the 
emitted power to increase superlinearly in the number of 
oscillators. The typical area over which this coherence 
can be realized is given by the pumping spot diameter. 
For example, a diameter of 60 jam would give an emit- 
ted power / ^ 11 jiiW. The emission intensity can be 
further enhanced by growing additional stacks of double 
QWs and placing the system inside a supplemental THz 
cavity. The latter would allow stimulated emission and 
a THz laser giving competitive characteristics compared 
to other solid-state THz sources. 

Conclusions. We have shown that nonlinear interac- 
tions in a dipolariton system give rise to multistability 
effects. In particular, for certain values of pumping pa- 
rameters the parametric instability between IX and DX 
modes occurs. This results in continuous oscillations of 
the spatially indirect exciton occupation number under 
CW pumping. The frequency of these oscillations is in 
the THz range and can be tuned by an applied electric 
field. Depending on the parameters, the resulting THz 



radiation by the array of classical dipoles represents a 
continuous single- mode or multi- mode THz laser, with 
power output and efficiency expected to be improved 
over existing solid state THz emitters. We have also 
shown rapid switching of the THz emission controlled 
by an additional short laser pulse. 
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